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The hadronic structure function of the photon F^(x, Q 2 ) is measured as a function 
of Bjorken x and of the photon virtuality Q 2 using deep- inelastic scattering data taken 
^ by the OPAL detector at LEP at e + e~ centre-of-mass energies from 183 to 209 GeV. 

Previous OPAL measurements of the x dependence of are extended to an average 
Q 2 of (Q 2 ) = 780 GeV 2 using data in the kinematic range 0.15 < x < 0.98. The Q 2 
evolution of F% is studied for 12.1 < (Q 2 ) < 780 GeV 2 using three ranges of x. As 
predicted by QCD, the data show positive scaling violations in with F^(Q 2 )/a = 
(0.08 ± 0.02 torn) + (0-13 ± 0.01 torn) ln Q 2 , wh ere Q 2 is in GeV 2 , for the central x 
region 0.10-0.60. Several parameterisations of F^ are in qualitative agreement with the 
measurements whereas the quark-parton model prediction fails to describe the data. 
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1 Introduction 



Much of the present knowledge of the structure of the photon has been obtained from mea- 
surements of the photon structure function F% in deep-inelastic electron-photonQ scattering 
at e + e~ colliders, see ]lj for a recent review. The large statistics and high electron energies 
of the full LEP2 programme permit the extension of the measurement of to higher values 
of (Q 2 ) than have been probed at LEP1. The photon structure function K7 is expected to 
increase only logarithmically with Q 2 p|. Therefore, the large range of Q 2 values accessible 
at LEP, which extends from about 1 GeV 2 to several thousand GeV 2 , makes it an ideal place 
to study the evolution. 

The measurement of F% in e + e~ interactions is based on the deep-inelastic electron- 
photon scattering reaction, e(k) "f(p) — > e(k') + hadrons, proceeding via the exchange of a 
virtual photon, j*(q), where the symbols in brackets denote the four-momentum vectors of 
the particles. The flux of quasi-real photons can be calculated using the equivalent photon 
approximation ||. The cross-section for deep inelastic electron-photon scattering is expressed 



as: 



d 2 o- C7 _ >cX 2ita 2 



(l + (1 - yf) F2(x, Q 2 ) - y 2 F2(x, Q 2 



dxdQ 2 xQ 4 

where Q 2 = —q 2 . The usual dimensionless variables of deep inelastic scattering, x and y, 
are defined as x = Q 2 /2(p ■ q) and y — (p • q)/{p • k), and a is the fine structure constant. 
The structure function F% is related to the charge-weighted sum of the parton densities of 
the photon (see e.g. |J). In the kinematic region of low values of y studied (y 2 1) the 
contribution of the term proportional to the longitudinal structure function F^(x, Q 2 ) is 
negligible 0. 

The analysis presented here is based on 632 pb _1 of data at e + e~ centre-of-mass energies 
^/s^ of 183 to 209 GeV, with a luminosity weighted average of ^/s^ = 197.1 GeV, taken 
by the OPAL experiment in the years 1997-2000. It extends the measurements of as 

1 For conciseness positrons are also referred to as electrons. 
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a function of x up to (Q 2 ) = 780 GeV 2 , and significantly improves on the precision of the 
measurement of the Q 2 evolution of F%. This analysis not only tests perturbative QCD but 
also measures F% at large Q 2 , a previously unexplored region in e + e~ collisions. This is 
approximately the region which has also been probed in jet production at HERA 0,[5|]. 

The paper is organised as follows. After the description of the OPAL detector in Section ^| 
the data selection is detailed in Section [| followed by the description of the Monte Carlo 
simulation and background estimates in Section |j. The results are presented in Section ^[ 
These comprise: the quality of the description of the observed hadronic final state by the 
Monte Carlo models, Section |5.1| ; the measurement of F% at high Q 2 , Section |5\2] ; and the 



measurement of the Q evolution of Section |5.3|. Conclusions are given in Section 



2 The OPAL detector 



A detailed description of the OPAL detector can be found in 0, and therefore only a brief 
account of the main features relevant to the present analysis will be given here. 

The central tracking system is located inside a solenoidal magnet which provides a uni- 
form axial magnetic field of 0.435 T along the beam axisQ The magnet is surrounded 
by a lead-glass electromagnetic calorimeter (ECAL) and a hadronic sampling calorimeter 
(HCAL). Outside the HCAL, the detector is surrounded by muon chambers. There are sim- 
ilar layers of detectors in the endcaps. The region around the beam pipe on both sides of 
the detector is covered by the forward calorimeters and the silicon-tungsten luminometers. 

Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector a precision vertex drift chamber, a large- volume jet chamber 
with 159 layers of axial anode wires, and a set of z chambers used to improve the measurement 
of the track coordinates along the beam direction. The transverse momenta p t of tracks 
with respect to the z direction of the detector are measured with a precision of cr Pt /p t = 

^/o.02 2 + (0.0015 • pt) 2 (p t in GeV) in the central region, 9 > 753 mrad. The jet chamber 
also provides energy loss, dE/dx, measurements which are used for particle identification. 

The ECAL covers the complete azimuthal range for polar angles that satisfy 9 > 200 mrad. 
The barrel section, which covers the range 9 > 609 mrad, consists of a cylindrical array of 
9440 lead-glass blocks with a depth of 24.6 radiation lengths. The endcap sections (EE) 
consist of 1132 lead-glass blocks with a depth of more than 22 radiation lengths, covering 
angles in the range 200 < 9 < 609 mrad. The electromagnetic energy resolution of the EE 
calorimeter is about 15%/\/E (E in GeV) at polar angles above 350 mrad, but deteriorates 
closer to the edge of the detector. 

2 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis upwards 
and the z axis in the direction of the electron beam. In this paper the polar angle 9 is defined with respect 
to the closest orientation of the z axis. 
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The forward calorimeters (FD) at each end of the OPAL detector consist of cylindrical 
lead-scintillator calorimeters with a depth of 24 radiation lengths divided azimuthally into 16 
segments. The electromagnetic energy resolution of the FD calorimeter is about 18%/y/E (E 
in GeV). The clear acceptance of the forward calorimeters covers the range 60 < 6 < 
140 mrad. Three planes of proportional tube chambers at 4 radiation lengths depth in the 
calorimeter measure the directions of showers with a precision of approximately 1 mrad. 

The silicon tungsten detectors (SW) || at each end of the OPAL detector lie in front of 
the forward calorimeters. Their clear acceptance covers a polar angular region between 33 
and 59 mrad. Each SW calorimeter consists of 19 layers of silicon detectors and 18 layers of 
tungsten, corresponding to a total of 22 radiation lengths. Each silicon layer consists of 16 
wedge-shaped silicon detectors. The electromagnetic energy resolution is about 25%/\/E (E 
in GeV). The radial position of electron showers in the SW calorimeter can be determined 
with a typical resolution of 0.06 mrad in the polar angle 9. 



3 Kinematics and data selection 



The interactions of two photons are classified according to the virtualities of the photons. For 
this analysis photons with a virtuality of less than 4.5 GeV 2 are called quasi-real photons, 
7, and the other photons are virtual photons, 7*. As a shorthand, events caused by the 
interactions of the three possible combinations are called 77, 7*7 and 7*7* events. 

To measure F^(x, Q 2 ), the distribution of 7*7 events in x and Q 2 is needed. These 
variables are related to the experimentally measurable quantities W, E tag and # tag by 

Q 2 = 2£ b £ tag (l-cos# tag ) (2) 

and 



Q 2 + W 2 + P 2 w 

where E h is the energy of the beam electrons, E tag and # tag are the energy and polar angle 
of the deeply inelastically scattered (or 'tagged') electron, W 2 is the invariant mass squared 
of the hadronic final state, and P 2 = —p 2 is the negative value of the virtuality of the quasi- 
real photon. The requirement that the electron associated with the quasi-real photon is not 
seen in the detector (anti-tag condition) ensures that P 2 <C Q 2 , so P 2 is neglected when 
calculating x from Equation [3]. The electron mass is neglected throughout. 

Three samples of 7*7 events are studied in this analysis, classified according to the subde- 
tector in which the scattered electron is observed. Electrons are measured using the SW, FD 
and EE detectors. Events are selected by applying cuts on the scattered electrons and on the 
hadronic final state. A scattered electron is selected by requiring E tag > 0.75/0.75/0.70 E h 
and polar angles 33.25/60/230 < # tag < 55/120/500 mrad for the SW/FD/EE samples. For 
the SW sample the energy cut effectively eliminates events originating from random coinci- 
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dences between off-moment umQ beam electrons faking a scattered electron and untagged 77 
events 0. For the EE sample special measures have to be taken to avoid fake electron candi- 
dates. To remove electron candidates originating from energetic electromagnetic calorimeter 
clusters stemming e.g. from hadronic final states in the reaction Z°/7* — > hadrons, an isola- 
tion cut is applied which requires that less than 3 GeV is deposited in a cone of 500 mrad 
half-angle around the electron candidate (electron isolation cut). 

To ensure that the virtuality of the quasi-real photon is small, the highest energy elec- 
tromagnetic cluster in the hemisphere opposite to the one containing the scattered electron 
must have an energy E a < 0.25 (the anti-tag condition). To reject background from 
deep-inelastic scattering events with leptonic final states, the number of tracks in the event 



passing quality cuts [|10j and originating from the hadronic final state, N tr ^, must be at 
least three/three/four for the SW/FD/EE samples, of which at least two tracks must not be 
identified as electrons, based on the energy-loss measurement in the jet chamber. The tracks 
and the calorimeter clusters are reconstructed using standard OPAL techniques |10| which 



avoid double counting the energy of particles that produce both tracks and clusters. The 
visible invariant mass W V i S of the hadronic system is calculated from tracks and calorimeter 
clusters, including contributions from energy measured in the SW and FD calorimeters. For 
the EE sample, because of the high probability that the scattered electron will shower in the 
dead material (ranging from 2 — 6 radiation lengths) in front of the EE calorimeter, energy 
deposits close to the electron are likely to belong to the electron. Therefore, for this sample, 
all tracks and clusters within a cone of 200 mrad half-angle about the direction of the elec- 
tron candidate are excluded from the calculation of W V i S . To remove the region dominated 
by resonance production and to reject the background from Z°/7* — > hadrons, the measured 
W vis is required to be in the range 2.5 < W vis < 60/60/50 GeV for the SW/FD/EE samples. 
The stronger cut on W vis applied to the EE sample reflects the fact that the background 
from Z°/7* — > hadrons is larger for this sample than for the other samples. 

The cuts applied to each sample are listed in Table |l|. The numbers of events in each 
sample passing the cuts are listed in Table |^, together with the numbers of signal events 
after subtracting the background contributions described below. Trigger efficiencies were 
evaluated from the data using sets of separate triggers, and were found to be larger than 
99% for events within the selection cuts. 



4 Monte Carlo simulation and background estimation 

Monte Carlo programs are used to simulate signal events and to provide background esti- 
mates. All Monte Carlo events are passed through the OPAL detector simulation [|ll[] and 
the same reconstruction and analysis chain as used for real events. 

The Monte Carlo generators used to simulate signal events are HERWIG 5.9+A: t (dyn) ||12|| , 

3 Off- momentum electrons originate from beam gas interactions far from the OPAL interaction region and 
are deflected into the detector by the focusing quadrupoles. 
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PHOJET 1.05 pi and the Vermaseren program [13]. The main reason for using a second 
Monte Carlo together with HERWIG is to have an additional model that contains different 
assumptions for modelling the hard scattering and the hadronisation process. HERWIG 
is a general purpose Monte Carlo program which includes deep inelastic electron-photon 
scattering. The HERWIG 5.9+k t (dyn) version uses a modified transverse momentum, k t , 
distribution for the quarks inside the photon for hadron-like events. The upper limit of the kt 
distribution is dynamically (dyn) adjusted according to the hardest scale in the event, which 
is of order Q 2 . This version was found to better describe the observed hadronic final states in 
three of the LEP experiments |15 than the original version HERWIG 5.9. In HERWIG the 



cluster model is used for the hadronisation process. PHOJET simulates hard interactions 
through perturbative QCD and soft interactions through Regge phenomenology, and the 



hadronisation is modelled by JETSET Hlg] . Since it is recommended by the authors to use 
PHOJET only for Q 2 values smaller than about 50 GeV 2 , the Vermaseren model is used for 
the EE sample. The Vermaseren program is based on the quark-parton model (QPM) and 
the quark masses assumed in the event generation are 0.325 GeV for u, d, s and 1.5 GeV for 
c quarks. For each Monte Carlo sample the generated integrated luminosity is at least 10 
times that of the data. 



The HERWIG and PHOJET samples were generated using the leading order GRV [17 
parameterisation of i^, taken from the PDFLIB library as the input structure func- 
tion. This version assumes massless charm quarks. Since PHOJET is not based on the 
cross-section formula for deep inelastic electron-photon scattering, the program always pro- 
duces the same x and Q 2 distributions independent of the input structure function. Therefore 
the x distribution of PHOJET was reweighted to match that from HERWIG, as described 
in [[HJ. This is not a strong limitation, because the main emphasis lies on the alternative 
hadronisation model. The result of the unfolding procedure is expected to be almost in- 
dependent of the actual underlying x distribution of the Monte Carlo sample used. The 
numbers of expected signal events from the HERWIG program are listed in Table [| 

For the SW and FD samples the dominant background comes from the reaction e + e _ — > 
e + e~T + T~ proceeding via the multiperipheral diagram |l[]. This was simulated using the 
Vermaseren program. In contrast, for the EE sample the dominant background stems from 
the reaction Z°/7* — * hadrons, which was simulated using PYTHIA f2"Uf . The next largest 
backgrounds are e + e~ — > e + e~T + r~ followed by non-multiperipheral four-fermion events 
with eeqq final states (denoted by 4-fermion eeqq), which were simulated with GRC4f [^1 



and Z°/7* — > t + t , which was simulated with the KK [^2 program. Because the aim is to 



measure the structure function of the quasi-real photon, events stemming from the interaction 
of two virtual photons with hadronic final states are also treated as background. For the 
SW and FD samples these were generated using PHOJET 1.10 with the virtualities of both 
photons restricted to be above 4.5 GeV 2 . For the EE sample they have been estimated using 
the Vermaseren program. The contribution to the background due to all other Standard 
Model processes was found to be negligible in all the samples. The numbers of events from 
the dominant background sources for each data sample are listed in Table |[ 
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5 Results 



5.1 Comparison of data and Monte Carlo 



Monte Carlo samples are used in an unfolding procedure to extract the differential cross- 
section da/dx and F%(x, Q 2 ) from the data. Therefore, apart from explicit effects due to 
the variation of the structure function, a good description of the data distributions by the 
Monte Carlo is needed both for electron variables, which are used to measure Q 2 , and for 
hadronic variables, which determine W 2 . The analysis of the SW sample closely follows that 



presented in ]0J but includes three times the data integrated luminosity The quality of the 
description of this sample is similar to that presented in |19j . The analysis of the FD and 
EE samples at LEP2 energies is new. Figures [I] and |2] show comparisons between data and 
Monte Carlo distributions for these two samples. The quantities shown are (a) E tag /E^, the 
energy of the scattered electron as a fraction of the energy of the beam electrons, (b) # tag , 
the polar angle of the scattered electron, (c) iV tr k, the number of tracks originating from the 
hadronic final state, and (d) W Y - m , the measured invariant mass of the hadronic final state. 
The FD sample, Figure [I], is compared to the Monte Carlo prediction of the HERWIG and 
PHOJET (without reweighting of the x distribution) signal events together with background 
estimates. The PHOJET sample has been normalised such that the predicted number of 
events for the SW sample is the same as that of HERWIG, so the PHOJET distributions 
only allow for a shape comparison. The HERWIG Monte Carlo model predicts slightly fewer 
events than are observed in the data and, in general, the shapes of the data distributions 
are better described by HERWIG than by PHOJET. The EE sample, Figure ||], is compared 
to the Monte Carlo prediction of the HERWIG and Vermaseren signal events together with 
background estimates. The data distributions of the energy and polar angle of the scattered 
electron are well described by the Monte Carlo predictions. For the variables related to the 
hadronic final state there are apparent differences in shape. 



The hadronic energy flow for the SW sample has been studied in |19||. On average, about 



5% of the energy is deposited in SW, and about 20-25% in FD and SW combined. The 
numbers for the FD and EE sample are even lower. It was verified that scaling the energy 
in the forward region has a small impact on the measured for x > 0.1. Consequently, in 
the present analysis this energy is not scaled. 

The quantity a; v i S , obtained from Q 2 and W v i S , is shown in Figure |3| for the three samples. 
It should be noted that for the SW sample, for x vis > 0.1, the HERWIG model using the GRV 
parameterisation of F% qualitatively follows the data, which means that the F% found from 
the data should be similar to the expectation from GRV. In contrast, for the FD sample for 
0.1 < x v i S < 0.7 the HERWIG prediction is systematically lower than the data. Due to the 
shortcoming of the PHOJET model discussed above, the description of the x v ; s distribution 
is unsatisfactory when using the PHOJET model without reweighting of the x distribution. 
For the EE sample the difference in shape of the W v - ls distribution is reflected in the observed 
difference between the data and both Monte Carlo models for the x vis distribution. 
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5.2 Measurement of dcr/dx and at high Q 

The differential cross-section da/dx and the structure function K7 are obtained from the 
data by unfolding the x vis distribution of the EE sample, after applying additional cuts on 
Q 2 . The main problem in measurements of at low x, i.e. x < 0.1, is the dependence 
of F2 on the Monte Carlo modelling, which enters when the unfolding process is used to 
relate the visible distributions to the underlying x distribution. This problem is less severe 
at medium to large values of x, in particular for the high Q 2 EE sample, where the hadronic 
final state has much more transverse momentum and as a consequence is better contained in 
the detector. Therefore the correlation between the measured invariant mass W V i S and the 
true W, e.g. as given by HERWIG, is much better at large x, so the results can be expected 
to have a smaller dependence on the Monte Carlo modelling of the hadronic final state. 

No attempt has been made in this analysis to access the region of x < 0.1, so using a 



one dimensional unfolding on a linear scale in x is appropriate, in contrast with |T9[. For 



this purpose the RUN program [^3|] has been used. Technically, RUN uses a set of Monte 



Carlo events which are based on an input F^(x, Q 2 ) and carry the information about the 
correlation of x v ; s and x. A continuous weight function is defined which depends only on 
x. This function is constructed from individual weight factors for each Monte Carlo event. 
These weight factors are obtained by fitting the x v i S distribution of the Monte Carlo sample 
to the measured a; vis distribution of the data, such that the reweighted Monte Carlo events 
describe the x v i S distribution of the data as well as possible. After the unfolding the two 
x vis distributions are consistent. The unfolded F^x.Q 2 ) from the data is then obtained by 
multiplying the input F%(x, Q 2 ) of the Monte Carlo with the weight function. For further 



details the reader is referred to M. It has been demonstrated in [IS] that this procedure is 



independent of the input structure function used in the Monte Carlo. 

Radiative corrections and the dependence of F^(x, Q 2 , P 2 ) on P 2 are treated as in the 
previous OPAL analysis [0J. The radiative corrections applied to the data have been esti- 
mated using the RADEG program |[24||. They are obtained for each bin in x and Q 2 using 



the SaSID |25|] prediction of F2- No correction for the effect of non-zero P 2 has been made, 
see Refs. [fj],[19| for further details. The average value of P 2 of the data samples as predicted 
by the HERWIG program is about 0.2 GeV 2 . Note however that HERWIG does not take 
into account the P 2 dependence of F2 . 

After subtraction of background, the EE sample has been unfolded using three bins in 
x spanning the range 0.15 — 0.98 and for 400 < Q 2 < 2350 GeV 2 . The central values are 
obtained using HERWIG as the input Monte Carlo model for the unfolding. Each data 
point is corrected for radiative effects as described above. Bin-centre corrections are also 
applied as given by the average of the GRSc [p6|, SaSID and WHIT1 |27| predictions for the 



correction from the average F% over the bin to the value of F^ at the nominal x position. 
The result for F% /a is shown in Figure f| and listed in Table |3] together with the correlation 
matrix. In each bin of x the result for da/dx is also listed. The da/dx values are corrected 
to the phase space given by the Q 2 range and y < 0.3. 
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Systematic errors are estimated by repeating the unfolding with one parameter varied at 
a time and determining the shift in the result. The systematic errors are combined by adding 
all individual contributions in quadrature separately for positive and negative contributions. 
The systematic effects considered for the EE sample are: 



1. Model dependence: 

The dependence on the Monte Carlo model used in the unfolding has been estimated 
by repeating the unfolding using the Vermaseren sample and taking the full difference 
as the systematic error, both for the positive and negative error. 

2. Variations of cuts: 

The composition of the selected events was varied by changing the cuts one at a 
time. The size of the variations reflect the resolution of the measured variables and 
the description of the data by the Monte Carlo models around the cut values. The 
variations are sufficiently small not to change the average Q 2 of the sample significantly. 
The variations made are listed in Table [l]. 

3. Unfolding parameters: 

The number of bins used for the measured variable can be different from the number 
used for the true variable. The standard result has 5 bins in the measured variable. 
This was in turn reduced to 4 and increased to 6 to estimate the systematic effects of 
the unfolding. 

4. Calibration of the tagging detector: 

The energy of the scattered electron in the Monte Carlo samples was conservatively 
scaled by ±1% p§. 



5. Measurement of the hadronic energy: 

The main uncertainty is in the calibration of the response of the electromagnetic 
calorimeter to hadronic energy for low energy particles in the hadronic final state. 
The absolute energy scale was varied by ±3% in the Monte Carlo samples. 

6. Background modelling: 

To quantify the uncertainty on the most dominant background, stemming from the 
reaction Z°/7* — > hadrons, the KK program along with cluster fragmentation from 
HERWIG has been used instead of PYTHIA with string fragmentation. 

7. Cone size for the W vis calculation: 

The size of the exclusion cone for the W v - ls calculation of 200 mrad half-angle about 
the direction of the scattered electron has been varied by ±30 mrad. 



The size of the contributions to the error from the individual sources is similar and no single 
source is dominant. When combining all error sources, the total estimated systematic error 
is of the same order as the statistical error. 

The measured F^/a, shown in Figure |] together with several theoretical calculations, 
exhibits a flat behaviour. The leading order parameterisations of from GRSc, SaSID and 
WHIT1, which all include a contribution from massive charm quarks, are described in detail 
in The contribution from bottom quarks is negligible. It can be seen that in this high 
Q 2 regime the differences between these predictions are moderate, particularly in the central 
x-region. All these predictions are compatible with the data to within about 20%, with the 
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WHIT1 parameterisation, which predicts the flattest behaviour, being closest to the data. 
The QPM curve, which models only the point-like component of F£, is calculated for four 
active flavours with masses of 0.325 GeV for light quarks and 1.5 GeV for charm quarks. 
This prediction shows a much steeper behaviour in x and is disfavoured by the data. 



5.3 Measurement of the Q 2 evolution of 

Following the study of the scaling violation of performed in |J0| the evolution of F^ 
with Q 2 has been measured for several x ranges using all three samples. Due to their large 
statistics, the SW and FD samples are further split into two bins of Q 2 (9-15 and 15- 
30 GeV 2 for SW and 30-50 and 50-150 GeV 2 for FD). The data are unfolded as a function 
of x separately in each bin of Q 2 and corrected for radiative effects. The results are shown 
in Figure |5] and listed in Table f|. The estimation of the systematic errors for the EE sample 
is described above. For the SW and FD samples the estimation of the systematic error 
mirrors the procedure for the EE sample, with some differences. The PHOJET program is 
used as a second Monte Carlo to determine the model dependence; the variations of cuts are 
given in Table [I]. For the unfolding parameters the standard number of bins, which was 8 
for the central values, has been varied by ±2. No systematics due to the electron isolation 
are needed for the SW and FD samples. For these samples, the largest contribution to 
the systematic error generally stems from the estimated model dependence. With only two 
models available that satisfactorily describe the data Hl9|1 , the estimated systematic error is 
small for those x, Q 2 regions where the two models happen to predict similar correlations 
between x and x v i S . To reduce fluctuations within the SW and FD samples, the systematic 
error from this source has been averaged for each region of x for the two Q 2 points within a 
given sample. 

The data in Figure |5](a) show positive scaling violations in i 7 "^ for the x ranges 0.10-0.25 
and 0.25-0.60. The QCD inspired parametrisations of F^ qualitatively follow the data, but 
do not perfectly account for them. For the SW sample the GRSc and SaSID predictions, 
which are almost indistinguishable, closely resemble the data, whereas at higher Q 2 , for the 
FD sample, WHIT1 comes closest. For the range 0.60-0.85, the data are compatible with 
the predicted scaling violations of the QCD inspired parametrisations. The QPM model 
generally gives a bad description of the data, especially at low x. 

To quantify the slope for medium values of x, where data are available at all values of 
Q 2 , the data are fitted using essentially the procedure from | 3~0fl . A linear function of the 



form a + b In Q 2 , where Q 2 is in GeV 2 , has been fitted to the data in the region 0.10-0.60. 
Within this range of x the parameters a and b are assumed to be independent of x. To 
obtain the central values of the two parameters, with their statistical errors and correlation, 
a fit was performed by the MINUIT [|31j program using the measured values of F^/a and 
their statistical errors as listed in Table (|. The fit was repeated for each of the systematic 
variations. The systematic errors of a and b are estimated as the quadratic sum of the 
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deviations of the two parameters from the central values. The result of the fit is 

F2(Q 2 )/a = (0.08 ± 0.02 + (0.13 ± 0.01 In Q 2 , 

with, for the central result, a correlation between the two parameters of —0.98 and a x 2 of 10 
for 3 degrees of freedom. No significant change of the result is observed if the fit is performed 
using the full error on each point. This new result compares to the previous OPAL value p0| 
of 

F2(Q 2 )/a = (0.16 ± 0.05±g;i2) + (0.10 ± 0.02 i™*) ^ Q 2 ■ 

These two determinations, based on independent data sets, are in agreement, and the errors 
on a and b have been significantly reduced. The data, together with the fit result, are shown 
in Figure 0(b). They are qualitatively described by the higher order GRV parametrisation 
(GRV HO). 



6 Conclusions 

The photon structure function K7 and the differential cross-section da/dx have been mea- 
sured using deep inelastic electron-photon scattering events recorded by the OPAL detector 
during the years 1997-2000 with an integrated luminosity of 632 pb _1 and an average e + e~ 
centre-of-mass energy of 197.1 GeV. 

The structure function F2 has been measured as a function of x in the range 0.15 < x < 
0.98 and at an average photon virtuality of (Q 2 ) = 780 GeV 2 , which represents the highest 
value measured so far. The Q 2 evolution of F^ has been studied for 12.1 < (Q 2 ) < 780 GeV 2 
using several ranges of x. The data exhibit positive scaling violations in i 7 ^ for the x ranges 
0.10-0.25 and 0.25-0.60. For the range 0.60-0.85, the data are compatible with the predicted 
scaling violations. The measured evolution of F^ ja as a function of Q 2 in the central region 
of x, 0.10-0.60, has been fitted with a linear function in InQ 2 , resulting in 

F2{Q 2 )/a = (0.08 ± 0.021°;°!) + (0-13 ± 0.01 l°;°i) ^ Q 2 , 
where Q 2 is in GeV 2 . 

Both for the measurement of F% at (Q 2 ) = 780 GeV 2 and for the investigation of the 
Q 2 evolution of F£ , the quark-parton model prediction is not in agreement with the data. 
It shows a much steeper rise than the data as a function of x for (Q 2 ) = 780 GeV 2 and 
also a different behaviour in the Q 2 evolution. In contrast, the leading order GRSc, SaSID 
and WHIT1 parameterisations and the higher order GRV parameterisation of i 7 ^ are much 
closer to the data. This means that the corresponding parton distribution functions of the 
photon are adequate to within about 20% at large values of x and at (Q 2 ) scales of about 
780 GeV 2 . 
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Cut \ Sample 


SW 


FD 


EE 


Et a „/Eh min 


0.75 (±0.05) 


0.70 (±0.05) 


ft™ min fmrad] 


33.25 (+2) 


60 (±3) 


230 (±5) 


^tae max [mrad] 


55 (-2) 


120 (-3) 


500 (±5) 


E a /E b max 


0.25 (±0.05) 


iVtrk min 


3(+l) 


4(±1) 




(2 non-electron tracks) 


W vis min [GeV] 


2.5 (±1) 


W v i s max [GeV] 


60 (±5) 


50 (±5) 


Electron isolation [GeV] 




3.0 (±0.5) 



Table 1: The selection cuts applied to each data sample, together with the variations applied 
(in brackets). See the text for explanation of the variables. 





SW 


FD 


EE 


data selected 


27819 


11874 


414 


data signal 


26071 ±167 


10652 ±110 


274 ±21 


Monte Carlo selected 


28308 ±51 


11211±32 


436 ±6 


HERWIG signal 


26560 ±49 


9989 ± 30 


296 ±5 


Backgrounds 


7*7 — > T + T~ 


1309.3 ±14.1 


845.5 ±11.3 


31.8±2.2 


7*7* — > hadrons 


321.3 ±4.7 


193.4 ±3.7 


5.0 ±0.3 


Z°/7* — > hadrons 


82.8 ±2.4 


124.6 ±3.1 


76.2 ±2.4 


Z°/ 7 * -> T+T- 


7.9 ±0.3 


10.5 ±0.4 


10.6 ±0.4 


4-fermion eeqq 


27.0 ±0.9 


48.2 ±1.1 


16.6±0.7 



Table 2: The numbers of selected events and signal events (selected events corrected for 
background) in the data compared to the signal predictions from the HERWIG program. 
The expected numbers of background events for the dominant sources according to Monte 
Carlo are also listed. The errors given are only statistical. 



16 



X 


X 


Fq/a 


da /dx 


radiative 


bin-centre 


range 


bin-centre 




[Pb] 


cor. [%] 


cor. [%] 


0.15-0.40 


0.275 


0.93 ±0.10lg;^ 


o.94 =b o.ii ±8:ii 


-8.8 


-0.4 


0.40 - 0.70 


0.550 


0.87 ± 0.10±g;?i 


0.79 ±0.09lg;?4 


-6.9 


0.1 


0.70-0.98 


0.840 


0.97 ±0.171^ 


0.62 iO.lllg;^ 


-5.4 


-5.1 



x range 


0.15-0.40 


0.40 - 0.70 


0.70-0.98 


0.15-0.40 


1.00 






0.40 - 0.70 


0.16 


1.00 




0.70 - 0.98 


-0.04 


-0.15 


1.00 



Table 3: Results for the EE sample for F^jo. as a function of x at (Q 2 ) of 780 GeV 2 , and for 
da/dx in the Q 2 range 400-2350 GeV 2 and y < 0.3. The first errors are statistical and the 
second systematic. The data were unfolded in bins defined by the x ranges and corrected for 
radiative effects. For a given bin the radiative correction is the difference of the radiative and 
non-radiative cross-sections as a percentage of the non-radiative cross-section. The structure 
function was corrected to the x values listed using the bin-centre corrections which are given 
as a percentage of the non-corrected F 2 7 . The statistical correlations between the bins for 
the central result are also given. 
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(a) 



x range 


Q 2 range 

1 Vjrt: V I 


(Q 2 ) 

\QpV 2 ] 

IVjrt: V J 


Fq/a 


da/dx 

[PDJ 


radiative 

rnr C%1 
LUI . [/OJ 


0.10-0.25 
0.25-0.60 


9-15 


12.1 


o.38 ±0.01 
0.43 ±o.oiig;g| 


83 ± 1 ± I 
47 ± 1 1 2 


-5.3 
-4.4 


n 1 n n 9^ 
0.25-0.60 


i ^ — 


1 Q Q 


u.oy ± u.ui _ 003 
0.49 ±0.01 ±8:82 


36 ± 1 1 J 


— 0.0 
-4.5 


n 1 n — o ?5 

0.25-0.60 
0.60-0.85 


30 — 50 


39 7 


47 +nni +0.02 

u.^ 1 n= u.ui _ 0.02 
0.63 ±0.01 ±8:8s 
0.65 ±0.03t8:8e 


01 7 4. n 4 + 0.9 

15.9 ± 0.4 1 8:f 
9.6 ± 0.4 1 8| 


— 5 9 
-4.8 
-3.8 


0.10-0.25 
0.25-0.60 
0.60-0.85 


50 - 150 


76.4 


0.55 ±0.01 18:8s 

0.68 ± 0.01 ±8:§1 
0.73 ±0.02 ±8:81 


18.8 ± 0.4 + 8:5 
13.8 ± 0.3 1 81 
9.1 ± 0.3 1 8i 


-6.5 
-5.2 
-4.1 


0.25-0.60 
0.60-0.85 


400 - 2350 


780 


0.94 ±0.09 toil 
0.83 ±0.11 ±8:29 


0.91 ± 0.09 ±8:?o 

0.71 ± 0.09 1 8is 


-7.6 
-6.0 



x range 


0.10-0.25 


0.25-0.60 


0.60-0.85 


0.10-0.25 


1 






0.25-0.60 


0.00/0.28/0.45/0.40/- 


1 




0.60-0.85 


-/-/-0.23/-0.19/- 


-/-/0.32/0. 27/0.17 


1 



(b) 



x range 


Q 2 range 


(Q 2 ) 


Fj/a 


da/dx 


radiative 




[GeV 2 ] 


[GeV 2 ] 




[Pb] 


cor. [%] 


0.10-0.60 


9-15 


12.1 


0.41 ± 0.01 +8:81 


57 ± 1 1 % 


-4.7 




15-30 


19.9 


o.46 ±o.oil8:81 


42 ± 1 1 \ 


-4.9 




30-50 


39.7 


o.58 ±o.oil8:8l 


17.7 ± 0.3 1 8:? 


-5.3 




50 - 150 


76.4 


o.64 ±o.oil8:8l 


15.3 ± 0.3 1 °:| 


-5.6 




400 - 2350 


780 


0.90 ±0.0918:1? 


0.86 ± 0.08 i°:l? 


-8.2 



Table 4: Results for the average /a in bins of x for several values of (Q 2 ) and for da/dx 
for several ranges of Q 2 . Section (a) corresponds to Figure |5](a), section (b) to Figure ||(b) 
The first errors are statistical and the second systematic. The data were unfolded in bins 
defined by the x and Q 2 ranges listed. The radiative corrections which have been applied are 
also listed. For a given bin the radiative correction is the difference of the radiative and non- 
radiative cross-sections as a percentage of the non-radiative cross-section. The statistical cor- 
relations between the bins for the central result for the (Q 2 ) values 12.1/19.9/39.7/76.4/780 
GeV 2 are also given. 
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Figure 1: Comparison of data distributions for the FD sample with Monte Carlo predictions. 
The open histograms are the sum of the signal prediction and the contributions of the major 
background sources, shown both for the HERWIG (full lines) and PHOJET (dotted lines) 
models. The Monte Carlo predictions are normalised to the data luminosity, except for 
PHOJET where the sample has been normalised such that the predicted number of events 
for the SW sample is the same as that of HERWIG. All selection cuts have been applied, 
except for any cut on the variable in the plot (indicated as dashed lines if within the region 
shown). The distributions are: (a) E tag / E^, the energy of the scattered electron as a fraction 
of the energy of the beam electrons; (b) ^tag, the polar angle of the scattered electron; (c) iV tr k, 
the number of tracks originating from the hadronic final state; and (d) W V [ S , the measured 
invariant mass of the hadronic final state. The errors given are only statistical. 
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Figure 2: Comparison of data distributions for the EE sample with the Monte Carlo predic- 
tions. See Figure [I] for details. 
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X vis 



Figure 3: The measured x vis distributions for the three samples (a) SW, (b) FD and (c) 
EE. The data are compared to Monte Carlo predictions containing signal and background 
contributions normalised to the data luminosity, except for PHOJET where the sample has 
been normalised such that the predicted number of events for the SW sample is the same as 
that of HERWIG. The errors given are only statistical. 
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Figure 4: The measured F^/at as a function of x for the EE sample. The data are unfolded for 
(Q 2 ) = 780 GeV 2 and compared to the leading order predictions from the GRSc (full line), 
SaSID (dotted line), WHIT1 (dashed line) and QPM (dot-dashed line) parameterisations of 
F2 /a. The inner error bars represent the statistical errors and the outer error bars represent 
statistical and systematic errors added in quadrature. The tick marks at the top of the figure 
represent the bin boundaries. 
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Figure 5: The evolution of /a as a function of Q 2 for several bins of x, (a) 0.10-0.25, 0.25- 
0.60 and 0.60-0.85 and (b) for the central region 0.10-0.60. The inner error bars represent the 
statistical errors and the outer error bars represent statistical and systematic errors added 
in quadrature. In (a) the data are compared to the predictions from the GRSc (full line), 
SaSID (dotted line), WHIT1 (dashed line), and QPM (dot-dashed line) parameterisations 
of F^/a. In (b) GRSc has been replaced by the higher order prediction from GRV and, in 
addition, the result of the fit is shown. 
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